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Dynamics of ultra-intense circularly polarized solitons under inhomogeneous 
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The dynamics of the ultra-intense circularly polarized solitons under inhomogeneous plasmas are examined. 
The interaction is modeled by the Maxwell and relativistic hydrodynamic equations and is solved with fully 
implicit energy-conserving numerical scheme. It is shown that a propagating weak soliton can be decreased and 
reflected by increasing plasma background, which is consistent with the existing studies based on hypothesis 
of weak density response. However it is found that ultra-intense soliton is well trapped and kept still when 
encountering increasing background. Probably, this founding can be applied for trapping and amplifying 
high-intensity laser-fields. 

PACS numbers: 52.38.Kd, 41.75.Jv, 52.35.Mw, 52.59.-f 
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The interaction between lasers and plasmas has at- 
tracted much attention since the invention of chirped 
pulsed amplification technology, 1,2 which makes ultra- 
intense laser with intensity (To > 10 18 W/cm 2 ) applica- 
ble. Interaction between such intense lasers and plasmas 
exhibits a rich variety of interesting nonlinear phenom- 
ena. One such phenomenon, namely the dynamics of 
electromagnetic solitons, has drawn much attention an- 
alytically, numerically and experimentally. 3-32 As nearly 
30% ~ 40% of the laser pulse energy can be transformed 
into solitons This fairly high efficiency of electromag- 
netic energy transformation indicates that solitary waves 
can play an important role in the development of the in- 
teraction between the laser pulses and plasmas. In fact, 
many theoretical problems are still open even for one- 
dimensional solitons, namely the properties of linearly 
polarized solitons^ and the soliton dynamics in the pres- 
ence of plasma inhomogeneity 2 ^— . In reality, the back- 
ground plasma density is always inhomogeneous, thus it 
is of great necessity to investigate the dynamics of soli- 
tons under inhomogeneous plasma background. 

Researches on the propagation of solitons under in- 
homogeneous background plasma have been conducted 
by various authors i^ 2 . - — The pioneering work in this 
field is on the theoretical studies of the Langmuir soli- 
ton through an inhomogeneous plasma. Their ana- 
lytical conclusions have demonstrated the acceleration 
(deceleration) of solitons in the presence of decreasing 
(increasing) plasma background. 22 However there is a 
obviously shortcoming in this work, the group veloc- 
ity of the soliton is treated as an independent param- 
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eter. In fact, the velocity, frequency and amplitude 
of the soliton are closely related, which is determined 
by an eigenvalue condition.— Later on, the propaga- 
tion of electromagnetic solitonic structures under inho- 
mogeneous plasma has been investigated i 23 ' 24 However, 
all these related researches have been limited to the 
hypothesis of weak density response.— The weak den- 
sity response theory can not deal with ultra-intense soli- 
tons accompanied with heavy plasma density cavitation. 
The dynamic behavior of solitons under rare background 
plasma density has been studied by 2D Particle-in-Cell 
simulations, 25 which shows that the solitons are acceler- 
ated toward the vacuum-plasma interface with an accel- 
eration proportional to the gradient of the plasma den- 
sity. Recently, fluid-Maxwell simulation has been car- 
ried out to systematically study the dynamic behavior 
of solitons under inhomogeneous plasma^ They have 
arrived at the conclusion— there is no qualitative differ- 
ence in the propagation features of large- amplitude solu- 
tions and the small- amplitude solutions of the nonlinear 
Schroedinger equation soliton variety, through the inho- 
mogeneous plasma background. They have found that 
both weak and strong solitons can be decreased and re- 
flected by the increasing plasma background. In their 
simulations, the propagating soliton is put in as an ini- 
tial condition, where the amplitude of the strong soliton 
is a ~ 0.9 (a = eE /m e cuJo) and that of weak soliton 
is a ~ 0.4. What about the ultra-intense soliton with 
amplitude a > 1.0, can its behavior still be the same? 

In this paper, the dynamic behavior of ultra-intense 
solitons under inhomogeneous plasma background has 
been conducted by fluid-Maxwell simulation which is 
solved with fully implicit energy-conserving numerical 
scheme. It is shown that the weak propagating soliton 
can be decreased and reflected by the increasing plasma 
background, which is consistent with the existing stud- 
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ies and can be described by the weak-density response 
theory. However the ultra-intense electromagnetic soli- 
ton can not be decreased and reflected by the increasing 
background, but is trapped and kept still in the plasma 
density cavitation. 

Compared with existing studies, we address the issue 
of propagation of electromagnetic soliton under inhomo- 
geneous background in a more general context, (a) The 
hypothesis of weak density response is dropped, and we 
can investigate the dynamics of ultra-intense solitons. (b) 
The background plasma is not restricted to be of rare 
density, instead, over-dense plasma background is con- 
sidered, (c) The soliton is self-generated by the interac- 
tion between laser and plasma as a more general process, 
which is not put in as an initial condition. 

The governing equations are the relativistic hydrody- 
namic fluid equations coupled with the full system of 
Maxwell equations. As the electron quiver velocity is 
much greater than the thermal velocity in relativistic 
intense laser light, we treat the plasmas as cold liq- 
uid. The normalized full system of equations are shown 
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where t = oj pe t, x — uj pe x/c, a ViZ — eA y ^ z /m e c 2 , 
n e ,i = n e j/n , p eA = p e>i /m e c, and 4> = ecj)/m e c 2 arc the 
normalized time, space, wave vector potential, electron 
and ion density and momentum (along the laser prop- 
agation direction x), and electrostatic potential of the 
space charge field, respectively. Furthermore, m is the 
ion/electron mass ratio, and u> pe and uq = m e uj pe / Aire 2 
are the plasma frequency and its related plasma density. 

To solve the equations, we use a fully implicit energy- 
conserving numerical scheme^^ The boundary condi- 
tions are a VtZ = 0, (f> = and n e ,i — at both sides of 
the simulation box. Our model is a global simulation, 
namely, at initial time, the incident laser and plasma are 
apart form each other. The laser impacts the plasma, 
and the evolution of vacuum-plasma interface is demon- 
strated in detail^ Unlike the existing studies that the 
propagating soliton is put in as an initial condition, in 
our simulation the soliton is self-generated through the 
interaction between laser and plasma. It is generally 
accepted that at vacuum region the fluid treatment of 



plasma is invalid, because the very definition of a plasma 
calls for an accumulation of macroscopic number of par- 
ticles. Numerical integrating of Eqs. (4) and (5) result 
in unphysical results of p e> i at the zero plasma density 
region. In our simulation, p e ^ is restricted to be zero at 
the vacuum region. As the zero plasma density leads to 
zero current density J x = 0, it has no effect on the solu- 
tion of coupled Maxwell Eq. (1), which makes the global 
fluid-Maxwell simulation applicable^^ 
The initially incident laser is of the form, 



a m cos 2 [(x — x p q)tt/2t] cos[cjo(£ — x)], 
a m cos 2 [(x — x p o)tt/2t} sin[wo(t — a;)], 



(8) 



where x p is the position of peak laser intensity, x p — r < 
x < x p + T, t = 5.0Ao (Ao = 2tt/u)q) is the laser duration, 
ujq = 0.88 is the laser frequency and a m = 0.71 is the 
maximum laser amplitude. The initial plasma density 
distribution is of the form, 

n e = n p sin 2 [(a; - x p i)ir/2(x p2 - x p i)], 

Xpi ^ X ^ Xp2) 

n e = n p + n p a sin 2 [(x - x p3 )ir/2(x pi - x p3 )], 

$p3 ^ % — «^p4 ; 

n e = n p (l + a) cos 2 [(x - x pi )Ti/2(x p5 - 3^4)], 

x p 4 < x < x p §. (9) 

To describe an increasing plasma background, we have 
n p = 1.0 (= 1.29n c ), a = 0.2, x p i = 0.0A o , x p3 = 4.0A 
and x P 4 — 8.OA0 The steepness factor is define as 
s = x p 2 — x p i , which plays an important role in the laser 
plasma interaction. As shown in Fig. HJa), the pene- 
tration ratio is closely related to s, the steepness of the 
vacuum-plasma interface. The penetration ratio is de- 
fined as {ei n — ebk)/£in, where 6j„ is the initial laser en- 
ergy and ebfc is the reflected light energy. We can see 
that a rather smooth vacuum-plasma interface allows for 
a strong penetration of laser pulse. According to this 
relation, we can control the amplitude of the soliton by 
adjusting the steepness of vacuum-plasma interface. For 
a rather steep vacuum-plasma interface, as shown in Fig. 
[TJb), most of the incident laser is reflected and the pen- 
etrated electromagnetic soliton is weak. However for 
a rather smooth vacuum-plasma interface, as shown in 
Fig- [He), the penetrated electromagnetic soliton is ultra- 
intense with amplitude a ~ 1.7. The behavior of ultra- 
intense soliton is beyond the prediction of weak-density 
response theory, and it is of great interest to numerically 
study its dynamics under inhomogeneous plasma back- 
ground. 

Figure [5] shows the space-time evolution of plasma den- 
sity (left column) and laser amplitude (right column) un- 
der inhomogeneous background. For the small-amplitude 
soliton with a ~ 0.5, as clearly demonstrated in Fig. [2] 
(al), the plasma density is almost undisturbed. Weak 
density response theory can well deal with the behav- 
ior of such small-amplitude soliton under inhomogeneous 
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FIG. 1. (color online) (a) The penetration ratio of the incident 
laser pulse vs. shape factor s. (b) and (c) The reflected 
and penetrated laser structure for different vacuum-plasma 
interface steepness, s = O.IAo and s = 0.4Ao, respectively. 
The red curve is for the normalized laser amplitude a = (a y + 
a 2 ,) 1 / 2 and black curve is for the plasma electron density. Here 
the simulation parameters are a y = a z = 0.7f, cjo = 0.88, 
r = 5.0Ao, x p o — — 5.02Ao and x p i = O.OAq. 



plasma background. The soliton undertakes decelera- 
tion and eventual reflection when encountering increas- 
ing plasma background, as shown in Fig. [2] (a2). As 
for the more intense soliton with a ~ 0.8, the plasma 
density is more or less distributed and a shallow cav- 
itation is formed moving with the propagating soliton, 
which is shown in Fig. [2] (bl). Strictly speaking, the 
weak density response treatment is no longer proper to 
describe the dynamic behavior of soliton with a ~ 0.8, 
however its behavior, deceleration and reflection, is al- 
most the same as the small-amplitude soliton. It seems 
that there is no qualitative difference in the propaga- 
tion features of large-amplitude solutions and the small- 
amplitude solutions of the nonlinear Schroedinger equa- 
tion soliton variety, through the inhomogeneous plasma 
background. However, Fig. [5] (c) and (d) indicate that 
large- amplitude solitons with a ~ 1.7 behave dramati- 
cally different from that of small-amplitude. The plasma 
density is strongly affected by the ponderomotive force 
of large-amplitude soliton, and a deep density cavitation 
is formed. Unlike the behavior of small-amplitude soli- 
ton, the large-amplitude soliton firstly experiences decel- 
eration, then is well trapped and kept still in the deep 
plasma cavitation, when encountering increasing plasma 
background . 

For a general wave packet with group velocity (3 prop- 



agating in the plasma, we have /3 = (1 



> 2 ) 



2U/2 



(1 — rtp/ne) 1 / 2 , 2 ^ where n p = 1.0. This relation demon- 
strates that the group velocity is closely related to the 
plasma background density. At n zero = 1/(1 — /3 2 ) = 
1/(1 — 0.17 2 ) = 1.03, the group velocity of soliton is de- 
creased to be zero. For small-amplitude soliton with al- 
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FIG. 2. (color online) The space-time evolution of plasma 
density and laser amplitude under increasing plasma back- 
ground, where the column (1) is for plasma density and 
column (2) for laser amplitude. Different cases of vacuum- 
plasma steepness (or soliton amplitude) with s = O.IAo, 
s = 0.2A , s = 0.4Ao and s = 0.8A (or a = 0.5, a = 0.8, 
a — 1.7 and a = 1.7) are shown in (a), (b), (c) and (d), respec- 
tively. Here the simulation parameters are a y — a z — 0.71, 
loq = 0.88, r = 5.0Ao, x p o = — 5.02Ao, x v \ = O.OAo, 
x p s — 4.0Ao, x P 4 = 8.OA0, n p = 1.0 and a = 0.2. 



most no density response, the group velocity of soliton 
is totally determined by the plasma background density 
profile. However for large-amplitude soliton, such as Fig. 
QJc), the density of the cavitation edge can be as high 
as 1.2, which is even higher than the increasing plasma 
background density. The high-amplitude soliton is sur- 
rounded by high plasma density located at its both front 
and back sides, which are both higher than n zero . In 
this case, the high-amplitude soliton can not be reflected, 
as the small-amplitude soliton does, by the increasing 
plasma background density. This high-amplitude and 
well trapped soliton can be the basis for trapping and 
amplifying high-intensity laser-fields between two solid 
density foils <2i 
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FIG. 3. (color online) (a) and (b) The space-time evolution of 
plasma density and laser amplitude under decreasing plasma 
background, (c) The back curve represents the incident laser, 
green curve represents the initial plasma density profile, and 
red curve for the penetrated laser pulse, (d) The spectral 
intensity of the incident laser and penetrated laser, where the 
black curve is for incident laser and red curve is for penetrated 
laser. Here the simulation parameters are a y — a z — 0.71, 
ujo = 0.88, t — 5.0Ao, x p o — — 5.02Aq, x p i = O.OAo, x P 3 = 
x P 4 = 4.0Ao, x p5 — 8.OA0, n p — 1.0 and s = O.8A0. 



Figure EJa) and (b) show the space-time evolution 
of plasma density and laser amplitude under decreasing 
plasma background. As expected, the soliton undertakes 
acceleration through the decreasing background. The 
electromagnetic soliton penetrates and enters the back- 
side vacuum of the plasma, where no second reflection 
occurs on the backside vacuum-plasma interface. As can 
be seen from Fig.[3jc), the amplitude and duration of the 
penetrated laser are stronger and shorter than that of the 
initially incident laser. Figure [3jd) shows the spectrum 
distributions of the incident and penetrated laser, which 
indicates that the penetrated laser is a good source of 
low-frequency burst of electromagnetic radiation. 

In summary, the dynamic behavior of ultra-intense 
soliton under inhomogeneous plasma background has 
been conducted by fluid-Maxwell simulation which is 
solved with fully implicit energy-conserving numerical 
scheme. The amplitude of propagating soliton can be 
well controlled by the adjusting the steepness of vacuum- 
plasma interface. It is shown that the small-amplitude 
soliton can be decreased and reflected by the increasing 
plasma background. However the large-amplitude soli- 
ton is well trapped and kept still encountering increas- 
ing plasma background. Probably, this founding can be 
applied for trapping and amplifying high-intensity laser- 
fields. The amplitude and duration of the penetrated 
laser pulse are much stronger and shorter than that of 
the initially incident laser. The penetrated laser is a good 
source of low-frequency burst of electromagnetic radia- 
tion. 
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